Phase fluctuations in anisotropic Bose condensates: from cigars to rings 
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We study the phase-fluctuating condensate regime of ultra-cold atoms trapped in a ring-shaped 
trap geometry, which has been realized in recent experiments. We first consider a simplified box 
geometry, in which we identify the conditions to create a state that is dominated by thermal phase- 
fluctuations, and then explore the experimental ring geometry. In both cases we demonstrate that 
the requirement for strong phase fluctuations can be expressed in terms of the total number of 
atoms and the geometric length scales of the trap only. For the ring-shaped trap we discuss the zero 
temperature limit in which a condensate is realized where the phase is fluctuating due to interactions 
and quantum fluctuations. We also address possible ways of detecting the phase fluctuating regime 
in ring condensates. 
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Since the laboratory realization of Bose-Einstcin con- 
densates in ultra-cold atomic systems, their properties 
in different dimensions and geometries and the effect of 
quantum and thermal fluctuations have been interesting 
questions. In Refs. [jj and Q, it was pointed out that 
for an elongated, 'cigar'-shaped condensate, a regime of 
strong thermal fluctuations can exist. Such a system is 
three-dimensional in the sense that its transverse dimen- 
sions are significantly larger than the healing length and 
the thermal de Broglie wavelength. However, the long- 
range behavior of the single particle (field-field) correla- 
tion function gi(r) = (^(O)^^)), with ip(r) being the 
single particle operator, can be dominated by phase fluc- 
tuations along the longitudinal direction. This leads to 
an exponential decay with some correlation length 1$, 
which is typical for a one-dimensional system at finite 
(non-zero) temperature. Because the many-body state 
combines these properties that are characteristic for 3D 
and ID systems, such a system can be considered to be of 
mixed dimensionality. When the correlation length 1$ is 
of the order of the system size or shorter, the condensate 
is said to be in the phase-fluctuating regime. In Refs. [3] 
and [4[ the 1D-3D cross-over was studied. In Ref. Q, a 
two-step condensation mechanism for the non-interacting 
system was pointed out. The phase-fluctuating regime of 
elongated condensates was first experimentally realized 
in Ref. jg|, and then studied in further detail in Refs. 
[7H10l| . We note that quantum fluctuations can lead to a 
phase fluctuating condensate even at zero temperature, 
as we discuss in Sect. Ill Bl However, here the decay 
of the correlation function is algebraic, and we find that 
the exponent is very small for the parameter regime of 
interest. The main focus of this article are thermal phase 
fluctuations. 

In this paper we study the properties of a ring-shaped 
phase-fluctuating condensate. One way to create such a 
system is by using both a quasi-2D optical dipole trap 
formed from a sheet of light (see Ref. Hi), and a 



Laguerre-Gauss (LG) beam perpendicular to the plane 
of the sheet [l^]. In such a setup, the atoms can form a 



condensate along the ring-shaped maximum of the com- 
bined intensities of the sheet trap and the LG beam. (A 
different approach was demonstrated in Ref. [13|.) Other 
realizations of toroidal traps were reported in Refs. [14j . 
[1 51 ] and [16j . Here we study under what circumstances 
the regime of thermal phase fluctuations can be reached 
for a condensate with a ring geometry. 

In Sect. QJ we first take the simplified case of a box 
geometry with periodic boundary conditions, with differ- 
ent box sizes along different axes, with one of them much 
larger than the others. In Sect. [II] we study the case of a 
ring-shaped condensate. In Sect. IHII we consider ways in 
which this state can be seen experimentally, and in Sect. 
HVlwc conclude. 



I. BOX GEOMETRY 

In this section we consider a condensate in a three di- 
mensional box with periodic boundary conditions, and 
dimensions X, Y, and Z. Our goal is to study the be- 
havior of the single-particle correlation function for the 
dimensional cross-over from three dimensions to one di- 
mension. This correlation function can be used to iden- 
tify the phase-fluctuating regime of the condensate. 

We consider the regimes Z « X, Y and Z ^> X, Y, 
using the calculational approach that is used in Refs. 1] 
and [2j. This Bogoliubov-de Gennes (BdG) approach is 
particularly useful for spatially inhomogeneous systems. 
However, we apply it to the box geometry as well, to allow 
for a clearer comparison between this case and the results 
for a ring geometry in Sect. [II] We first calculate the 
BdG modes, and then consider the thermal occupation 
of these modes. 

In a qualitative sense, the box geometry with periodic 
boundary conditions displays the main features of a Bose- 
Einstcin condensate in elongated geometries. Because it 
is calculationally simpler, it will serve as the main system 
to understand the properties of fluctuating condensates. 
In Sect. [II] we derive the analogous properties for the spa- 



tially inhomogeneous case, and compare them to results 
of this section. 

We consider an atomic ensemble of bosons of density 
n = N/V, where N is the total particle number, and 
the volume V is XYZ. These atoms are assumed to be 
weakly interacting via a contact interaction with strength 
g = 4nh 2 a/M, where a is the scattering length, and M 
is the atomic mass. We then apply a BdG approach to 
describe the condensate and the excitations. The con- 
densate is described by the well-known Gross-Pitaevskii 
(GP) equation, Ref. [lj: 
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2M 



ff|$o| *0 = M*0 
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where ^o is the condensate wave function, and /j, is the 
chemical potential. Within the Thomas-Fermi (TF) ap- 
proximation, see Ref. [19], ^o is given by 



*o = y/rio, 



(2) 



where n = N /V is the condensate density, and -/V 
is the total number of atoms in the condensate. For the 
weak interactions considered here, the quantum depletion 
is small and Nq ss N. The Bogoliubov-de Gennes (BdG) 
equations |18| are given by 

/ h 2 v 2 \ 

V ~2M~) Uu + 9\^o\ (2u v ~ v ") = (H + E v )u u (3) 
/ h 2 V 2 \ 

{-^f)v»+g\%o\ 2 (2v v -u v ) = (n-E„)v„.(4) 

The fields u v and v v describe an excitation of energy E v . 
We note that for a phase-fluctuating condensate, this for- 
malism should be interpreted as an expansion in the den- 
sity around a finite value, and in the phase to quadratic 
order, see Ref. [2J|. We introduce the rescaled quantities: 
x = x/X, y = y/Y, z = z/Z, and C 2 = h 2 /(2Mfj,L 2 ) 
where L is defined as L = (XYZ) 1 / 3 . Our objective is 
to generate a consistent expansion in £. Written in these 
rescaled quantities, the GP equation ([T]) is given by 
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The reduced density n is given by n — (7I \J>o| 2 /a*- We 
note that the TF solution is simply the zeroth order so- 
lution in £, i.e. uq = 1. The rescaled BdG equations 
are 



fluctuations of excitation v. In terms of these fields, the 
BdG equations §j\ and © arc 



-C 2 V 2 /+,„ = 2Ce„/_ 
< 2 V 2 /-,„ + 2/_^ = 2Ce r ,/ 4 



(9) 
(10) 



The solutions can be written as /± jm = C± tm W(m). 
Here we replaced the generic index v by the vector index 
m that appears in the solutions. The prefactors have to 
be related as C_ im = £ e m(7+ !m , due to the normaliza- 
tion conditions on u m and v m , see Refs. [1J and 0. The 
functions W(m) are plane waves, W(m) — exp(27rimf), 
where m = (m\ 1 m2, W3), and mi,m,2 and 7713 are inte- 
gers, f is defined as f = (x,y,z). Eliminating /_. m from 
Eqs. © and (I10[) . we can obtain the low-momentum 
limit of energy-moinentum dispersion relationship 
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Written in terms of the original parameters, this gives 
E m = hc\k\, where the phonon velocity is c = \/gn /M, 
and k is defined as k = (2nmi/X 1 2irm2/Y,2'Km3/Z). 
The normalization of the BdG modes is given by Cj_ m = 
2/j,/(E m V). In summary, the solutions of the BdG equa- 
tions are 
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Within the weak interaction expansion used here, we 
have [i = griQ. 

We now use these modes to calculate the correlation 
function of the phase. The phase field <j)(r) is related 
to /+, m through 4>(r) = l/(2y/n5) J2 m /+, m aL + h.c, 
where a^ is the bosonic creation operator of the /+. m 
mode. The /_. m modes, on the other hand, are related 
to the density fluctuations of the system. Because the 
modes /+, m behave as l/y/E m , whereas /-, m scale as 
y/E m , one can expect that at long distances (small k) 
the phase of the system can be a fluctuating quantity, 
whereas the density fluctuations will be suppressed. 

The correlation function of the phase is given by 

(50(r) 2 ) = ^-f^AU2-2cos(27rmr)), (14) 

m VEm 

where Scf>(r) = 0(0) — <fi(r). We assume that the Bogoli- 
ubov modes are occupied according to a Bose-Einstein 
distribution of temperature T. Thus, the distribution 
iV m , which is the occupation of a state with energy E m , 
is given by 



( 2 V 2 u„ + n (2u p 
n (2v l/ 
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u„) = (l-2Ce v )v u , (8) 



where 2(e u — E v /fi. We now define the fields f± iV = 
u v ±v„. These fields are related to the phase and density 



N n 
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(15) 



where the 1/2 accounts for the vacuum fluctuations. For 
small energies and at finite temperature, N m has a diver- 
gent behavior, which is the origin of the phase-fluctuating 



regime. It behaves as 



Nr, 



k B T 

Em 



(16) 



This divergent term dominates over the constant term 
1/2, which we discuss in Sect. Ill Bl We find that it 
generates an algebraically decaying correlation function, 
with an exponent that is very small for the parameter 
regime of interest. In particular, the term does not con- 
tribute to the exponential decay, and will be ignored in 
this section. 



A. Correlation function in ID and 3D 

The correlation function for the elongated geometry 
shows 3D behavior on short length scales r < X,Y, and 
ID behavior for longer scales r > X,Y. In this section 
we discuss these limits. 

We first calculate the 3D limit by replacing the sums 
in Eq. ([14]) with integrals. The thermal contribution of 
the correlation function is 



W>(r) 2 ) = 
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where r — |r|. For distances r large compared to the heal- 
ing length £ = h/(2Mgno) 1 ^ 2 and the thermal de Broglie 
wavelength A = h/(2nMk B T) 1 / 2 , this approaches 



(8<Kt) 



2ir 2 hc [ 



"K 2 (k B T) 2 TTk B T 
3h 2 c 2 her 
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(see e.g. Ref. [25J]). This expression shows the 1/r decay 
of the phase correlations of a 3D BEC to a non-zero value, 
related to the thermal depletion at T < /j,. The thermal 
depletion is proportional to £/(nA 4 ). Equation (fTS]) is 
valid in the elongated geometry for £, A <C r < X, Y. 

As we approach the ID limit, for r ^> X,Y, the modes 
with mi = rri2 = produce the dominant contribution 
to the correlation function, whereas the excited trans- 
verse modes provide corrections, which we discuss be- 
low. We write the correlation function as (S(f)(r) 2 ) = 
(S(f>(r) 2 )o + (Scj)(r) 2 )i, where (<50(r) 2 )o corresponds to the 
modes with mi = 7772 = 0, and (5(f>(r) 2 )i to the remain- 
ing contributions in Eq. (fT4|) . For (Scj)(r) 2 )o, with the 
approximation in Eq. (|16|) . we find 



W(r) 2 



gk B T Z 
XY Jhcj- 
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(19) 



In this expression, the linear term contains the length 
scale that we are interested in, whereas the quadratic 
term ensures that the derivative of the function at z = 
1/2 is smooth. One might expect an expression that is 
purely given by trigonometric functions of z, however, 
the approximation in Eq. (|16[) leads to the expression 
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FIG. 1: (Color online) Thermal contribution to the 
correlation function (?i(r) = (i/' t (0)V'(r)). We plot 
ln((^r (0)^(r))/no), which equals — (6</>(r) )/2. We use a box 
geometry with Z = 100/^m and X = Y = 5/jm, and a temper- 
ature T - lOOnK. The total atom number is N = 20000, and 
M is the sodium mass M = 23u, with u being the atomic mass 
unit. We assume that the scattering length is a — 2.75nm. 
The upper diagram shows the full correlation function, the 
lower diagram the regime z < X, Y . 



above. We note that if we take the limit Z — > 00, the 
linear term stays fixed, because it is only a function of z, 
whereas the quadratic term vanishes, because it scales as 
1/Z. 

At intermediate distances X, Y <C r <C Z (where the 
periodic boundary conditions are not apparent), the sin- 
gle particle correlation function thus approaches 



(^((Mr)) ~ 77 exp(-Z|z|/^), 



(20) 



where we used the phase-density representation of the 
boson operator ^(r) = y/n(r) exp(i0(r)) to connect the 
single particle correlation function to the phase correla- 
tion function. The phase correlation length Z^ is given 
by 



2XY{hc) 2 _ 2XYh 2 n 
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We note that the interaction has cancelled in this expres- 
sion (see e.g. Refs. [J] and j27|). This can be understood 
by writing the Hamiltonian in the phase-density repre- 
sentation. As discussed in Ref. [2J|, the Hamiltonian of 
a weakly interacting Bose gas is given to second order 
by a kinetic term ~ no(Vc/>) 2 /2M and a potential term 
of the form gSn 2 , where Sn are the density fluctuations 
around the constant value no- In the momentum repre- 
sentation, the spectrum of this Hamiltonian can be seen 
to be equivalent to the Bogoliubov spectrum of Eqs. ^ 
and (U]). The contribution to the partition function made 
by a given mode of energy Huik consists of independent 
phase and density terms, for a temperature T > tuuk/ks. 
In this limit, the phase correlations are entirely due to 
the kinetic term, leading to the expression in Eq. (|2T|) . 
independent of the interaction strength. 

Next, we estimate the effect of the higher modes, for 
the case X — Y. We consider the thermal contribution to 
the spatially independent part of the correlation function. 
We replace the 7773 summation in Eq. (|14[) by an integral 
over k z , and expand the dispersion relation appearing in 
the exponent in Eq. (fT5j) to second order in k z . We can 
then evaluate the k z integral as 
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Next we approximate the sum J2m u m 2 jto h V El=i 2? 
and (m 2 + m 2 ,) 1 ' 2 by m. That gives 
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where Li_i is the standard poly logarithmic function 
with exponent —1/2. 

For temperatures large enough to populate transverse 
modes, i.e. k B T 3> 2irhc/ X, but smaller than or compa- 
rable to /j,, Eq. ([2"B")) becomes 
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2 3 / 2 tt(>V) 3 ' 
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We therefore approximately recover the 3D thermal de- 
pletion, which we found before in Eq. (TTg|). For k B T <C 

2irhc/X, (S(j>(r) 2 )i is exponentially suppressed 



WW 2 ). 



2g (k B T\^/ 2 



hcXKXhc 



exp 



2ttHc \ 
Xk B TJ' 



(25) 



a) 








x^ 


T [nK] 1 






1 / 




0.75 








T 






/ 


-*2^ 


c 


0.5 




/ 
/ 
/ 


^^y' 


__,-' 


0.25 


fa " 






III 






1 


2 3 , 


4xl0 13 






n. [cm 3 ] 




b) 2 






L J 




T[nK] 

1.5 






111^ 


,--" T . 


1 






s ^^^--^^ 


^ff~ 


0.5 




^^ 




1 

■ 



N, 



4xl0 5 



FIG. 2: (Color online) Comparison of T c (black, continuous 
lines) and T^ (red, dashed lines), for different geometric aspect 
ratios. In a) we show these temperatures as a function of 
density, for the sodium mass M = 23u, and Z — 100/im, and 
for X = Y — 10, 7.5, 5/im (I - III). In b) we show the same 
quantities as a function of the total atom number No. T c is 
purely a function of the density, whereas T$ is a function of 
the total number, for a fixed longitudinal dimension Z. 



In Fig. [T] we show the full thermal contribution to 
the correlation function, Eq. (|14p. for box dimensions 
Z = 100/im and X = Y = 5/im. The temperature is 
T = lOOnK, and the atomic mass is M = 23u, with u 
being the atomic mass unit. For z < X, Y the correla- 
tion function 51 (r) closely resembles the 3D result (see 
Eq. (jTTJ)), i.e. it falls off on a short scale to a finite value, 
which describes the thermal depletion of the condensate. 
For sufficiently large z ^> X,Y, the correlation function 
<7i(r) is approximately given by gi(z), and behaves like 
the correlation function for a ID system. The thermally 
activated modes where 7711,7712 = 0, 7773 7^ generate an 
exponentially decaying g\{z). In this ID- like case, the 
correlation function is essentially only a function of the 
coordinate z, and independent of x and y. The effect 
of the transverse modes is to create a constant thermal 
depletion, as in Eq. (|2"4"]l . which is of the order of the 3D 
result. So the cross-over from 3D to ID is given by the 
geometric shape, and the regimes are z < X, Y for 3D 
behavior, and z > X, Y for ID behavior (see also Ref. 



B. Phase-fluctuating condensate 

In the previous section we found that for an elongated 
system the long range behavior of the single-particle cor- 
relation function 51 (r) shows exponential decay with a 
correlation length 1$. We can then define a tempera- 
ture above which the system is in the phase-fluctuating 
regime, based on the requirement that this length scale is 
comparable to the system size, i.e. l<f, ~ Z. This defines 
the following temperature 



knTit, = 



2h 2 N 
MZ 2 



(26) 



If this temperature lies below the BEC temperature T c , 
a phase-fluctuating condensate may exist. We estimate 
that temperature by using the result for the homoge- 
neous, ideal Bose gas, which is (e.g. Ref. [23j |) 



kuT r 



2ix 



2/3 



C(3/2) 2 / 3 M 
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where C, signifies the Riemann zeta function. To esti- 
mate the regime in which one would observe a phase- 
fluctuating condensate, we consider the ratio of T c and 
Trk, 



TJT^ 



1 Z 4 / 3 



C(3/2)V3 N i/axW 



(28) 



which we wrote in terms of the total atom number iVo, 
and we chose X = Y . For a system of fixed dimensions 
X and Z , no phase- fluctuating regime exists if the total 
number of atoms fulfills 



-7T- 3 7 4 7 4 

N„ > _JE ^«45— 

~ C(3/2) 2 X 4 ' X 4 ' 



(29) 



demonstrating how sensitively the presence or absence 
depends on the geometric ratio Z/X. In Fig. [5] we show 
both the condensate temperature and T^ as a function of 
density, for different geometric dimensions. The regime 
above T^ but below T c is the phase-fluctuating conden- 
sate regime. We can again see that the magnitude of this 
regime is very sensitive to the ratio Z/X. 



II. RING GEOMETRY 

We now turn to the kind of ring-shaped geometry in 
which we are interested. As mentioned in Sect. HI we 
will find that qualitatively most features of the phase- 
fluctuating condensate in a ring can be understood in 
terms of the simplified box geometry discussed previ- 
ously. We note that the elongated box with periodic 
boundary condition is very similar to a ring. In addition, 
the cigar-shaped condensates discussed in Refs. [l| and 
2] qualitatively share similar features in the center of the 
condensate. In this section, we derive the results analo- 
gous to those of the previous section but for a trapped 



(inhomogeneous) ring system and compare them to re- 
sults of Sect. HI 

In the treatment here, the trapping potential is ap- 
proximated by the harmonic form 



V 



-Mu 2 z z 2 + -Mu%p-R) 2 . 



(30) 



Here, z and p are spatial cylindrical coordinates as shown 
in Fig. [3l and ui z and uj p the oscillation frequencies along 
the spatial directions, and R is the mean ring radius. We 
apply the same formalism as in Sect. |U The GP equation 
of the system is given by 
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the BdG equations are 
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= (p + E v )u v (32) 
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= (p-E v )v v . (33) 
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The Thomas-Fermi (TF) solution of Eq. J3IJ) for z 2 
p 2 < 1 is 



#0 = V /n max(l-Z 2 -p 2 )- 



(34) 



Here we have rescaled the spatial coordinates z — z/l z 
and p = (p — R)/l p by the TF radii, given by l z = 
(2/i/Mw 2 ) 1 / 2 and l p = (2/i/Mw 2 ) 1 / 2 . This rescaling is 
in analogy to rescaling the spatial coordinates of the box 
geometry by the dimensions of the box. n max is the max- 
imum value of the density of the condensate. We write 
the GP equation as 

- 77 2 V 2 * + z 2 ^ + p 2 * + ^0*0 = *o (35) 
where we introduced the operator 
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Id Id 2 



uj 2 V dp 2 p dp 



u^d 2 
uj 2 dz 



a 8fa ,.(36) 



where 9 is the azimuthal angle, as shown in Fig. [31 and 
p = p/lp. The frequency uj = ^/uj z uj p , and rj 2 = lp |28l] . 

We now introduce a rescaled energy e v = E v jhw and 
a rescaled density no, defined as no = |\E'o| 2 /'ramax. With 
these parameters, the BdG equations are 

- rfV 2 u v + (z 2 + p 2 )u v + (2u v —v v )no 

= (l + 2r)e u )u u (37) 
-n 2 V 2 v u + (z 2 + p 2 )v u + (2v u - u v )ho 

= (1 - 2r)e u )v v . (38) 

These equations are analogous to Eqs. (J7J and ([5]), but 
here the equations contain a spatially dependent poten- 
tial and condensate density tt,q = n o(z,p). 




p = 'pp 



FIG. 3: (Color online) The toroidal shape of the BEC is given 
by the radius R, and the Thomas-Fermi radii l z in z direction, 
and l p in radial direction. We use the rescaled z variable 
z = z/l z , and the two rescaled radial variables p = p/l p and 
p — (p — R)/l p , as indicated in the figure. 



As before, we introduce the fields /±.„ = u v ± v v . In 
terms of these, the BdG equations are 
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(39) 



r/ 2 V 2 + 2 - 2z 2 - 2p 2 
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f-,u 



= 2r,e v f+, v (40) 

We write the solutions in the TF approximation as 

f+,u = c + Ai-z 2 -p 2 ) 1/2 w v (e) (4i) 

f-,u = C^^(\-z 2 -p 2 Y 1 l 2 W v {B) (42) 

with the normalization condition C_.„ = r)e u C+^. The 
functions W v (9) are simply W v = exp(im#), with m be- 
ing an integer. We replace the index v with m from now 
on. For these solutions we assumed that the radius R 
is much larger than l„, so that we can approximate the 
term {l/p 2 )d 2 in Eq. flU by (l 2 p /R 2 )d 2 . 

This set of modes is analogous to the modes with mi — 
rn 2 = in the previous section. In analogy to Eq. (|lip . 
the dispersion relation can be calculated from 



(1 - 5 2 - f 



>(■ 



V 2 ^o 



U 
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where we need to average over the density profile, see 
Refs. [lj and [20|. This gives E m = hc\k\, with c = 
\/gn ma , x /2M and k — mj R. These modes and disper- 
sion relation were also discussed within a hydrodynamic 
approach in Ref. [2fJ. The prefactor C+, m is given by 
C\ m = l/(2ir 2 e m 7iRl z l p ), so /±, m are 
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exp(im9). (45) 



These expressions are analogous to Eqs. (|T2"j) and (TT31 
The 'Thomas-Fermi' volume V in Eqs. (fT5|) and |T 



corresponds to 2ir 2 Rl z l p , which is the volume of the torus. 
The chemical potential p, is replaced by gn (p, z). 

Next, we calculate the correlation function of the phase 



mof 



E 



gN„ 



2ir 2 Rl z l p E m 



(2-2cos(m0)) (46) 



where ./V m is given in Eq. (|15[) . In analogy with 5(f>(r) in 
the previous section, we have defined 6(f)(0) = </>(0)- 
Using Eq. (|16p . we find 



WW 2 ) 



gksT R 

TT 2 l z l p (He)'- 



■{-9 2 /2 + n\9\), (47) 



in analogy to Eq. 19. So for intermediate dis- 
tances l p ,l z < i?0 < i? we have (^(9)^(9')) ~ 
exp(-i?|A(9|/^), with 



^r6 — 



'zi-p' 



fi 2 Wn 



Mk B T 



ixMRkeT 1 



(48) 



in analogy to Eq. 20. Here we used the relationship 
between the density maximum n max and the total atom 
number No, Nq — n max ir 2 l z l p R. Again we find that the 
interaction does not affect 1&, We note that this expres- 
sion is of the same form as the expression for the box 
geometry, Eq. (|2"Tj) . up to numerical prefactors. 



A. Phase-fluctuating regime 

As in Sect. HI the phase-fluctuating regime is reached 
when the phase correlation length is of the order of the 
system size. From 1$ ss ttR we can define the temperature 



kWTs = 



h 2 N 

ir 2 MR 2 ' 



(49) 



We calculate the condensation temperature for a cylin- 
der of length 27rf? and with harmonic confinement with 
frequencies uj z and w p 



knT r 



Ur(a 



3N h 3 u z ujp \ 2 /5 



Ar(a)C(a)(2m) 1 / 2 R 



(50) 



where C( a ) is again the Riemann function at a — 5/2, 
and r(a) is the usual Gamma function. 

As discussed in Ref. [2l|, at this critical temperature 
the density at the potential minimum reaches the critical 
value of the homogeneous system, within a semi-classical 
approximation. 

The ratio T c jT§ can be written as 



T A V 



3tt 5 



2/5 1 / R 2 \4/5 



V472r(a)C(a)/ 7V 3/5 Va P 



) i=wr(zz) < 51 > 



where a z = yjh/mu} z and a p — ^h/mu p are the oscilla- 
tor lengths of the confining potential. For the total atom 




25000 



FIG. 4: (Color online) T$ (red, dashed line) as a function of 
the total atom number, for M = 23m and R — 40^m, and 
T c (black, continuous lines) for u z = uj p — 2tv x 200, 2n x 
100, 2tt x 50Hz (I -III). 



to be of the order of the healing length lh 

.22]. 
The scaling exponent K is given by 



K 



HN, 



3/4 



rV2 



2 1 /4i?3/4 M 3/4 £,1/4(^^)1/4- 



see 



(55) 



For typical parameters, K is very large. For example, 



for R = 40^m, M = 23u, 



2.75nm, u> z — oj„ — 



2tt x 1000Hz, and N = 20000, we have K = 242.4. Since 
the algebraically decaying correlation function reaches, 
say, the value 1/e at a distance r c e 2K , which vastly ex- 
ceeds the size of the system, the decay of the correlation 
function across the system is very small. So although the 
system in principle forms a quasi-condensate, for typical 
parameters it will behave essentially like a true conden- 
sate. 



number this leads to the inequality 
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(52) 
(53) 



for which there is no phase-fluctuating condensate for 
any temperature. This conclusion is valid if the estimate 
of the condensate temperature in Eq. (J50I) is accurate, 
and under the weak interaction conditions of the BdG 
approach. 

In Fig. U]we show T^ and T c as a function of the total 
atom number iVo, and for different confining frequencies 
w 2 and uip. For a given Nq and trap radius R, T$ is fixed. 
The condensation temperature, however, is affected by 
the transverse confinement through uj z and lu p . By in- 
creasing their value, a large regime of a phase-fluctuating 
condensate can be created. On the other hand, for suf- 
ficiently weak transverse confinement, one can avoid the 
phase fluctuating regime for all temperatures. 



B. Zero temperature 

While one might expect that in 3D a condensate would 
have perfect phase coherence at zero temperature, we 
find that the phase-fluctuating character of an elongated 
system persists even at zero temperature, as for a truly 
ID system [29j . At zero temperature, the single-particle 
correlation function of our ring shaped condensate decays 
algebraically: 



(ROf 



<^(nMr 2 )) ~ (i^-) 



-1/4.R- 



(54) 



where ri and r 2 both are at the minimum of the trapping 
potential, with |r! | = |r 2 | = R, and 9 is the angle between 



them. Equation (jM)) is valid for 



■zi ^pi 



r c < \BB\ < R, 



III. EXPERIMENTAL SIGNATURES 

In this section we consider possible experimental sig- 
natures of the phase-fluctuating regime. The defining 
feature of this regime is the loss of phase coherence over 
the ring-shaped condensate. One way of testing this 
would be to use the interference technique demonstrated 
in Ref. [13. 

In Ref. [6fl the signature of the phase-fluctuating regime 
in time-of-flight images was studied. In Ref. [7[ an inter- 
ferometric approach measuring the second-order corre- 
lation function was reported. In Ref. [9] Bragg spec- 
troscopy was used to measure the momentum distri- 
bution in the axial direction, which was calculated in 
Ref. [8J . Matter wave interferometry was used in Ref. [10| 
to measure the phase coherence length. 

A consequence of the ring geometry is that during a 
time-of-flight (TOF) expansion of the phase-fluctuating 
condensate, the contrast of matter-wave interference 
should be reduced compared to a true condensate. In 
particular, due to the ring-shaped geometry of the system 
one can expect that the constructive interference that can 
emerge at the center of the ring during expansion will be 
reduced. 

We now discuss the interference of atoms released from 
a ring trap. For this purpose we use a Green's func- 
tion approach which conveniently relates the density after 
time-of-flight to the two-field correlation function prior 
to release. Assuming freely expanding atoms, the time- 
dependent bosonic single particle annihilation operator 
is given by 



■0(x,t) = / dxiw(x — Xi, t)-0(xi). 



(56) 



where r c is a short-range cut-off which should be chosen 



Here, rp(pc,t) is the bosonic single particle annihilation 
operator at time t, and i/'( x i) is f ne initial single parti- 
cle operator. The Green's function of free propagation 
w(x, t) is given by 

/ M \ 3 / 2 
™(x,i) = (3-^) exp(zAfx 2 /2ftt). (57) 




FIG. 5: (Color online) Expectation value of the density, in- 
tegrated over z direction, after time-of-flight expansion of an 
atomic ensemble of No = 30000 atoms of mass M = 23«, 
in a ring with radius R — 40/xm, and ui z = 2% x 1000Hz 
and uj p = 2-k x 200Hz. a)-c) show the TOF interference for 
T = 20nK, at times t — 10, 20, 30ms, d)-f) show the case of 
T = 170nK. 



The expectation of the density at time t is therefore 
n(x, t) = 

dxidx 2 w*(xi,i)w;(x2,t)(-0 t (xi)V'(x2)) (58) 



with Xi = x — xi and x 2 = x — x 2 . We are interested in 
the density integrated over the z direction, and we switch 
to cylindrical coordinates: 



n(p,i) = / dzn(p, z,i). 



(59) 



Using Eq. 
n(p,t) 



and integrating over z, we find 
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\2ttMJ M, 

exp(iM(p(pi cos Q\ — p 2 cos # 2 )) /fit 

-iM(p\ - pl)/2ht)y/n(pi,zi)n(p 2 ,zi) 
exp(-.R/^(|A0| - (A9) 2 /2n)). (60) 



Here we assumed that the single particle correlation func- 
tion is of the form 



(4>HP1, Zl,0l)*[>(P2, Z 2 ,02)) 



= y/n(p 1 ,z 1 )n(p 2 ,z 1 )exp — \^~ ~ \ A9 \) ( 61 ) 

with A9 = 9\ — 92- Here we have ignored the condensate 
depletion due to thermal activation of transverse modes. 
We perform the integral of Eq. (|M|) using shifted coor- 
dinates p\ and f>2, defined by p\ = R+pi and p 2 = R+P2- 
Ignoring quadratic terms in pi i2 in the exponent we ob- 
tain 



n(p,t) 



dz 1 d9 1 d9 2 A*(9 1 )A(9 2 ) 

exp(iMpR(cos9i — cos9 2 )/frt) 
exp(-R/lt(\A0\ - {A9) 2 /2n)) (62) 



where A{9) is given by 
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(63) 



We evaluate Eq. (|6"2")l numerically, for R — 40nm, 
7V = 30000, uj z = 2tt x 1000Hz, and uj p = 2tt x 200Hz. 
We assume M to be the sodium mass, M — 23u, where 
u is the atomic mass unit, and a — 2.75nm. We note 
again that in this example the interactions are taken into 
account for the initial state, whereas the time-of-flight 
expansion is ballistic. The initial velocity distribution is 
given by the TF wavefunction in the transverse direction 
and by the thermal distribution of the phononic modes 
along the longitudinal direction. 

In Fig. [5] (a) - (c) we show the case of T = 20nK, 
at times t — 10, 20, 30ms, in Fig. [S] (d) - (f) the case of 
T = 170nK for the same times. The latter is in the phase- 
fluctuating regime and the contrast of the expectation 
value of the density is visibly reduced. In particular, the 
density at the center varies with time as 



n(0,t) 



F(Z ) / dz{\ - z 2 /ll) 



J 2 (Ml p 1/I - z 2 /l 2 z R/M) 



(64) 



where 



F(h) 



si I 

-^ exp(-^/2/ )er/,( v / ^/2/ ) (65) 



is a measure of phase coherence around the ring. 

For small l^/R, Fil^) behaves as F(l^) « 21$/ R. For 
large values of l$/R, it saturates to F(^) = 27r. This 



behavior of the overall scale of the interference directly 
reflects the phase coherence of the BEC. The more co- 
herent the system is, the higher the interference contrast 
of the time-of-flight image. 

At short times n(0,£) shows a linear behavior in time, 
when time-averaged over the fast oscillatory behavior: 



n(0,i) 



At long times n(0, t) falls off as 
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htl. 
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(66) 



(67) 



By comparing Eqs. (|6"6")) and (|S7| we can see that the time 
at which the center density is maximum scales as t max ~ 
R/(H/(l p M)). The radius R is the distance traveled, and 
h/(l p M) is the scale of the velocity distribution in radial 
direction. 

We note that the expansion of a repulsively interact- 
ing gas should be faster than that of a non-interacting 
one, because the potential energy of the initial state is 
transformed into kinetic energy. The main observation, 
however, that a phase-coherent sample would show con- 
structive interference at the center of the ring, whereas 
an incoherent sample would not, should still apply. How- 
ever, the increased atomic velocity would lead to smaller 



fringe spacings making optical resolution of the narrow 
central peak more difficult. For sufficiently small in- 
teractions and for high enough optical resolution, this 
method could nevertheless constitute a method to detect 
the phase-fluctuating regime. 



IV. CONCLUSIONS 

We have studied the phase-fluctuating regime of con- 
densates in a ring-shaped trap. We first considered a 
simplified box geometry in Sect. HI and calculated the 
properties of the phase correlation function. We then 
turned to a realistic ring geometry in Sect. [II] The main 
results are the estimates of the T^ in Eq. (|49|) . and of T c 
in Eq. (J5T))) . With these expressions, one can determine 
whether the phase-fluctuating regime can be reached, as 
shown in Fig. 2J We present a simple condition, depend- 
ing only on atom number and trap aspect ratio, that 
guarantees avoidance of the phase fluctuating regime for 
all temperature below T c . We also discussed observing 
the regime of phase fluctuations in time-of-flight images. 
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